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A simulation method based on the RG blocking is shown to yield statistical errors smaller than that of the 
crude MC using absolute values of the original measures. The new method is particularly suitable to apply to the 
sign problem of indefinite or complex-valued measures. We demonstrate the many advantages of this method in 
the simulation of 2D Ising model with complex-valued temperature. 



1. THE SIGN PROBLEM 

In order to evaluate a multi-dimensional inte- 
gral of the partition function 



Z = fdV 



(1) 



using Monte Carlo (MC) one can sample the 
points in the integration domain with a non- 
uniform distribution, p. This sampling gives the 
following estimate for the integral: 



Z « (f/p)±y/S/N, (2) 

where N is the number of points sampled, p > 
and is normalised; 

1 N 

v/p) = jrYtf&yp^- 

i=l 

and 



S = j \f/ P -Z\ 2 pdV, 

« (f 2 /p 2 ) (f/p) 2 ■ 



(3) 



The best choice of p is the one that min- 
imises the standard deviation squared S. This 
can be found by variational method leading to 
the crude average-sign MC weight 



(4) 



Pcrudc = 1/1/ J \f\dV, 

giving the optimal 



■udc 



\f\dV 



fdV 



(5) 



The sign problem jij arises when Zj J \f\dV is 
vanishingly small: then unless a huge number of 
configurations are MC sampled, the large statis- 
tical fluctuations of the partition function render 
the measurement meaningless. 

2. IMPROVED SIMULATION METHOD 

One way of smoothing out the sign problem is 
to do part of the integral analytically, and the 
remainder using MC. The analytical summation 
is not just directly over a subset of the dynam- 
ical variables; in general it can be a renormali- 
sation group (RG) blocking where coarse-grained 
variables are introduced. This does yield certain 
improvement over the crude MC in general. 

Let P{V, V} be the normalised RG weight re- 
lating the original variables V to the blocked vari- 
ables V' 



P{V',V} > 0, 
J P{V',V}dV = 1. 



Inserting this unity resolution into the integral (Q) 
/ = JdV J dV'P{V',V}f, 



dV'g{V). 



(6) 



Thus, an MC estimator is only needed for the 
remaining integration over V' in (^|). As with 
the crude method of the last section, variational 



2 



minimisation for S of (|^), with g in place of /, 
leads to the improved MC 



Pimprovcd 



\g\dV . 



(7) 



Firstly, the improved weight sampling yields a 
partition function of magnitude not less than that 
sampled by the crude weight: 



(( Si g n ))impr 



= z 



\g\dV 
z 



J\J P{V',V}fdV\dV r ' 



> Z/J \f\dV . 

= l« si g n »crudcl 



(8) 



Secondly, it is also not difficult to see that the 
statistical fluctuations associated with improved 
MC is not more than that of the crude MC, 



S; 



improved 5crude 
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Pimprovcd 



d y-J 


f 




Pcrudc 



dV, 



WV ) (/ lfldV ) ' 



<0, 



(9) 



Thus the RG blocking always reduces the sta- 
tistical fluctuations of an observable measurement 
by reducing the magnitude of v3?/|((sign))|. 

Note that the special case of equality in (||,|9|) 
occurs iff there was no sig n problem to begin 
with. How much improvement one can get out of 
the new MC weight depends on the details of the 
RG blocking and on the original measure /. 

3. COMPLEX 2D ISING MODEL 

The Hamiltonian for the Ising model on a 
square lattice is 



H 



{nnl) 



(10) 



Here we allow j and h to take on complex values 
in general. For the finite lattice, periodic bound- 
ary conditions are used. 




Figure 1. Phase diagram in the complex j plane. 



The phase boundaries for the complex temper- 
ature 2D Ising model with h = are depicted in 
Figure 1 §. 

For the improved method, we adopt a simple 
RG blocking over the odd sites, labeled o. That 
is, the analytic summation is done over the config- 
uration space spanned by the even o sites; while 
MC is used to evaluate the sum over the remain- 
ing lattice of the odd • sites. In general, with 
finite-range interactions between the spins, one 
can always subdivide the lattice into sublattices, 
on each of which the spins are (nearly) indepen- 
dent for the partial sum to be carried out. 

Summing over the spins s Q , 



Z = 



JJ 2 cosh [js+ + h] (11) 



where s+ is the sum of the spins of the four odd 
neighbours surrounding each even site. The im- 
proved MC weight is then the absolute value of 
the summand on the right hand side of the last 
expression for Z. 

The quantities measured are magnetisation, 



3 



0.2 0.3 



Figure 2. | << sign >> | vs Re(j) Im(j) = 0, 
h = 0. Filled boxes are the improved data. 




M , and susceptibility, x- These can be expressed 
in terms of the first and second derivatives of Z 
respectively, evaluated at h = 0. In all the simu- 
lations, square two-dimensional lattices of various 
sizes with periodic boundary conditions are used. 
The heat-bath algorithm is used. Two additional 
benefits arise from the improved method. The 
first is that the number of sites to be visited is 
halved. While the expressions to be calculated at 
the remaining sites turn out to be far more com- 
plicated, the use of table look-up means that eval- 
uating them need not be computationally more 
expensive. The second benefit is that the number 
of sweeps required to decorrelate data points is 
reduced. 

As a test of the improved method, it is com- 
pared to the crude one along the path OX in 
Figure 1. Comparison with series-expansion data 
are plotted in Figure 3, in which the discs are the 
absolute-value circles of complex- valued statisti- 
cal errors for simulated results. 

4. CONCLUDING REMARKS 

The sign problem is partially eleviated in the 
improved method || because of some partial 
phase cancellation among the original indefinite 



Figure 3. Improved \ vs Re(J); expansion data 
shown as line. 



or complex-valued measure after an exact RG 
transformation. A particular RG blocking is cho- 
sen for our illustrative example of the 2D Ising 
model with complex-valued measure. But other 
choices of RG blocking are feasible and how ef- 
fective they are depends on the physics of the 
problems. 
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